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ABSTRACT 

A metric is introduced on the two dimensional space of parameters describing the Ising 
model on a Bethe lattice of co-ordination number q. The geometry associated with this 
metric is analysed and it is shown that the Gaussian curvature diverges at the critical 
point. For the special case q = 2 the curvature reduces to an already known result for 
the one dimensional Ising model. The Gaussian curvature is also calculated for a general 
ferro-magnet near its critical point, generalising a previous result for t > 0. The general 
expression near a critical point is compared with the specific case of the Bethe lattice and 
a subtlety, associated with the fact that the specific heat exponent for the Bethe lattice 
vanishes, is resolved. 
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§1 Introduction 

The concept of defining a distance function, or metric, on the space of states of 
a thermodynamic system or statistical distribution has recently been investigated by a 
number of authors (for a review from the statistical point of view of see [1] , from a 
thermodynamical point of view see [2] ) . The idea has also been extended to field theories 
[3] , [4] , [5] , [6] . For two dimensional field theories a metric was used to powerful effect in 
the proof of the c-theorem [3] . The purpose of this paper is to investigate further the role 
of a metric in thermodynamics and statistical mechanics by examining in detail a class of 
exactly soluble models — Ising models on a Bethe lattice of arbitrary co-ordination number, 
q. Near the critical point, the analysis should not depend in detail on the underlying model 
and it is shown how the specific results for the Bethe lattice fit in with a general analysis 
for a ferro-magnet near the critical point. 

A metric may be defined on the two dimensional space of states, which can be pa- 
rameterised by the spin-spin coupling J and the external field H, or by K = J/T and the 
magnetisation per site M. Near the critical point the curvature can be expressed purely in 
terms of critical exponents and the scaling function, provided the specific heat exponent 
< a < 1, and it will be shown how the general result for the curvature relates to that of 
the Bethe lattice, for which a = 0. This generalises the result of [2] for the curvature at, 
and just above, the critical temperature with vanishing magnetisation. 

In the specific example of the Bethe lattice, with co-ordination number q, it will be 
shown that the Gaussian curvature is positive and tends to q/2 for large temperature, but 
diverges at the critical point K c = jr = \ In (^2). These results generalise those of the 
one dimensional Ising model with periodic boundary conditions derived analytically in [7] , 
and are significant because they are exact results for a class of models which exhibits a 
phase transition. 

§2 The metric 

The motivation for a metric comes from probability theory and the theory of large 
deviations [8] . For example in thermodynamics [2] two states are considered to be close to 
one another if the probability of a fluctuation between them is large — if the probability is 
unity the distance between the states should be zero, if the probability is zero the distance 
between the states should be infinite. This is satisfied if the distance is a multiple of a 
positive power of minus the logarithm of the probability. The metric is thus related to the 
entropy. 

Let g a denote the intensive variables (temperature, pressure, external magnetic field, 
etc.) and <p a the corresponding extensive variables. The probability of a fluctuation from 
4> a to 4>' a , for given values of g a corresponding to a most probable value </> a , is given by a 
Gaussian distribution for 5(f) a = 4>' a — (j) a small, 

P ((/>', g) = Ne~ S{ ^' )+ ^~^^} HaHb+ 
2 



(1) 



where Af is a normalisation factor and S the thermodynamic entropy (the Einstein sum- 
mation convention is used throughout). 



Taking — In P we see that a natural metric is given by 

dl " = -ir^r d ^ d ^ ( 2 ) 

0(paO(pb 

or, in terms of the Legendre transform to W = — In Z 

W{g) = -S^) + g a $ a (3) 

with g a = 

* - -W^ b (4) 

We shall work on a lattice with iV sites and set w = = jp^, the free energy per site 
divided by the temperature (we use units in which Boltzmann's constant k = 1). It will 
then be convenient to rescale the line element (4) by the constant factor 1/N and use 

dl 2 = -J^-Sg a dg b . (5) 

QgaQ g b y V > 

The reason for this rescaling is in order to ensure that the metric makes sense in the 
thermodynamic limit. If the line element were not rescaled all distances would be infinite 
in this limit, since the probability of a fluctuation between two states would always be zero. 
The curvature associated with the line element (5) is finite in general, except possibly at 
critical points (this will be seen explicitly for the Ising model on a Bethe lattice in the 
following). Thus the curvature associated with the line element (4), which is 1/N times 
that of (5), would always vanish in the thermodynamic limit, except at critical points 
where it would be ill defined. This procedure of rescaling quantities by powers of N in 
order to study their behaviour under fluctuations is standard practice in the theory of large 
deviations. Thus the interesting geometry and non-zero curvature associated with (5) is 
due to thermal fluctuations. 

Equation (5) is not co-variant under general co-ordinate transformations, but that is 
because it is written in a special co-ordinate system. If 4> a = \^N < $ a > is the expectation 
value of the rescaled extensive quantity \^N^ a: then a co- variant form of the line element 
is 

dl 2 = {($ a $ b ) - ^ a )(^ b )}dg a dg\ (6) 

which reduces to (5) when the partition function is the sum over states of the exponential 
of a linear combination of \fN^ a with co-efficients g a , 

Z(9)= £ e-v^. (7) 
states 
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The form of equation (5) is only preserved under linear transformations of g a , not under 
general co-ordinate transformations. However equation (6), which reduces to (5) when the 
partition function is of the form (7), is general co-ordinate co- variant. In what follows it 
will be convenient to use the special co-ordinate system (7) in which the metric takes the 
form 

d 2 w 

QgaQgl 



For systems with two variables, e.g. K = ^ and h = ^ with J a spin-spin coupling 
and H an external magnetic field, it is convenient to perform the Legendre transform from 
(K, h) to (K, M) with M = — Changing to these variables, using 

dg c dg^ 

° a ' b ' = dg^dg^ ^ ^ 

in equation (8), one finds that the metric in the co-ordinate system {g a ) = (K,M) is 
diagonal, 

Ga* b i=(- T KK ° ), (10) 

\ U 1 MM J 

where 

T(K,M) = w + hM and T KK := ^ T M M := j^. (11) 
Equations (10) and (11) provide the basis of the calculations in the following two sections. 
§3 Curvature of the Bethe lattice 

The Ising model on a Bethe lattice, of co-ordination number q, has partition function 



Z = exp K Oi<jj + h Oi 



(12) 



where means a sum over all spin configurations on the lattice and Yl ^ s a sum over 

nearest neighbour pairs, K = J/T is the spin-spin coupling and h = H/T an external 
magnetic field. A good reference for a description of the properties of this model is [9] , 
whose notation is used here. For sites deep inside the lattice the free energy per site, /, is 
given by 

f/T = - X -qK - X -q ln(l - z 2 ) + ± ln[^ 2 + 1 - z(x + x~ x )] + 1 (q - 2) ln(x + x" 1 - 2z) (13) 

where z = e~ 2K and x is defined implicitly by 

_ e~ 2K + e - 2h xi- 1 
x ~ 1 + e -2h-2K xq -i ( 14 > 
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with < x < oo. This model has a phase transition at T c given by K c = J/T c = 
\ ln[q/q — 2]. The magnetisation per site is 



M = 



e 2h - x q 
e 2h + xi 



(15) 



which ranges from -1 for x — > oo to + 1 for a; = 0. A derivation of equations (13) — (15) 
starting from (12) can be found in [9]. Explicitly, equations (14) and (15) give 



x(M, z) 



y/1 - (1 - z 2 )M 2 + zM _ 1 

(1 + M) ~ x(-M, z) 



(16) 



The metric (10) is diagonal if the variables K and M are used, so consider the Legendre 
transform T = f/T + Mh 



F(K,M) = -q\n 



2z< 



{1-z 2 ) 



s-z \ 1 

M , 

+ Y ln 



1 — M 2 \ oM , 
, + V- In 

4 

1 + M 
1-M 



s + zM 
1 + M 



(17) 

where z = e 2K and s = y/l — (1 — 2 2 )M 2 . Equation (16) ensures that r(iiT, M) = 
T(K, —M). The derivation of equation (17) from equation (13) can be performed by using 
equations (14) and (15) to deduce 



dx 



dM 
dx 



dK 



x(x + x 1 — 2z) 2 
K ~~ ~ 2{2 - z(x + x- 1 )} 

2zx(x — x~ x ) 
m {2-zix + x- 1 )}- 



(18) 



(19) 



As described in the introduction, the metric adopted here involves the second derivatives 
of r 

-Fkk 

fjiM 



G 



(20) 



One finds, 



MM 



(y - izf 



Tkk = 8q 



(l + z< 



+ 



+ 



q 



zy) {zy - 2) 

z 



{1-z 2 ) 2 (y-2z)(zy-2) 



(21) 



where y = x + x 1 = 2 ^ 1 zM - • The Gaussian curvature can be determined 

from the metric (20) using (21). The algebra is tedious, but the final result is remarkably 
simple 

(1 - z 2 ) {2(1 + z 2 - yz) - q(3z 2 - 2zy + 1)} 
2 {(zy - 2) + q(l + z 2 - yz)} 2 ' 
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This can be written explicitly as a function of M and K, with z = e 2K and s = 

y/l-(l-Z 2 )M 2 , 

n( K, M) = V {(«-2)( fe -.)+4, } 
V ; 2 (s-z){(q-2)(^-s) + 2z} 2 V ; 

For T — > oo (if = 0) this reduces to 

K(0,M) = q/2 (24) 

and for M 2 = 1 

TZ(K,l)=q/2. (25) 
The curvature is plotted in figure 1 for g = 3, higher q give qualitatively similar pictures. 

Near the critical point, T = T c (l + t) with M << 1 and t << 1, z is of the form 
z = z c {l + 2K c t + o(t 2 )} with K c = \ ln(^) and z c = One then finds 

n „ 1 (g- 1 ) 2 1 f26) 

^2^-2) { i ln( _^_ ) , + (2 _i )M2} - 

In particular, 7?. diverges to +oo at the critical point (M = 0,t — > 0). 

The line of spontaneous magnetisation for t < is given by the smallest solution of 

see reference [9], and 1Z is finite along this line, except at the critical point M = t = 0. 
For small negative t, equation (27) gives the equation for the spontaneous magnetisation 

M 2 = (-t) J q \ - \J -*-) +0 (t 2 ). (28) 



2(q-l) \q-2 

Together with equation (26) this determines the Gaussian curvature along the line of 
spontaneous magnetisation as 

K=i4p^ ! s + ott\ (29) 

which is finite except at the critical point t = 0. 

The curvature diverges along the line on which the magnetic susceptibility diverges 
(the pseudo-spinodal curve [10] ). This line lies in the unstable region and does not coincide 
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with the line of spontaneous magnetisation (except when q — > oo) — though they touch 
at the critical point. These two curves are plotted in figure 2 for the case q = 3. 



These results are in complete agreement with the general analysis in [2]. Equation 
(23) reduces to the known result for one dimensional Ising model where q = 2 and M = 

sinhV\/sinh 2 h + e~* K , [7]. 

§4 Curvature Near the Critical Point 

Near the critical point the "effective action" for a general ferro-magnet 

T(t,M) = £ + Mh (30) 

is assumed to have the scaling form, independent of the microscopic model, 

T = \t\ 2 ~ a F(x) + analytic part, (31) 

where x = t/M 1 ^ is a scaling variable and F(x) a scaling function, a and f3 are standard 
critical exponents (see for example [11] ). In the variables t and M, the metric is diagonal 

/ o \ 

G = f X )• (*»> 
v u am ' 

Thus (—Gu) is the specific heat and Gmm the inverse magnetic susceptibility. 
In terms of the scaling function, F(x). 

G tt = -|t|- Q {(2 - a)(l - ct)F + (3 - 2a)F + F} 

U|2-a-2/3 . (33) 

G M m= " p2 x 2 ^(PF + F) 

where F = xjj^- etc. Stability requires Gu > and Gmm > 0, which puts conditions on 
the function F, 

(2 - «)(1 - a)F + (3 - 2a)F + F < 

x 2 P([3F + F) > 0. (34) 

In order that the internal energy U = — ^ be finite at t = 0, we must have a < 1, thus 
(2 — a) (1 — a) > 0. For x > eliminating F from (34) requires 

(3 - 2Q '-' 3) F<F. (35) 



(2-a)(l-a) 
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The calculation of the Gaussian curvature from equations (33) is straightforward but te- 
dious. The final result can be succinctly expressed by introducing a "dual" scaling function 
F{x) = x 2 ~ a F{x). The result is 



11 = 



(a + (3-1) 



(1 F 2 ) (2- a- p)(FF - F 2 )x 



2-a 



{T - T) {(3F + F) 



■ (36) 



(37) 



This extends the result of [2] to the case of non-zero magnetisation, and so is valid for 
t < as well as t > 0, (for comparison note that the result of reference [2] is quoted in 
terms of a different function, ^-Y(M/t 13 ) = F(t/Me)). The derivation of equation (36) 
assumes that the \t\~ a term in G tt (equation (33) ) dominates over the analytic part, i.e. 
it assumes that a > 0, so it is only valid for < a < 1. This means that equation (36) 
cannot be applied directly to the Bethe lattice at the critical point, for which a = 0. 

The scaling function, F(x) in equation (31), for the Bethe lattice can be obtained by 
Taylor expanding T(K, M) in equation (17) about M = and extracting the co-efficient 
oft 2 . First expand T(K,M), 

+ "6T\ 24 "^ (1 " z)3(322 + 9z + 8) 

+ o(M 8 ). 

Now let t = y c , so that z = e 1 + t where K c = ± ln^^^j, giving 

r = D + Et + t 2 F(x) + o(t 3 ) (38) 
where the scaling function has the form (also derived in reference [9]) 

F(x)=a+- + 4, (39) 
x x z 

and we have defined the following constants 



2 U(g-l) 



(40) 
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a = -\ j^j^ {2(3g - 2){q - 1) + g( g - 2)K C } 

6=^-2)K c (41) 

(remember a; = t/M 1 ^). The subtlety here is that, when a is not strictly positive, (32) 
is not quite correct because it assumes that the non-analytic term dominates over the 
analytic part, which is not true when a = 0. In general one has 

(dK\ 2 „ dfdKdY\ d 2 K 8Y 

G tt = hr g kk = -—(——) + 



V dt J dt\ dt OK J dt 2 dK , . 

\ / 42 

d 2 K(dK\~^ 



-~ Ttt + ywvdT) 

For a > T tt dominates over the second term on the R.H.S., but for the Bethe lattice 
a = 0, so the second term must be retained. Taylor expanding K = ^ gives 

G tt = -r tt -2r t + o(*). (43) 

Using (40), (39) and (41) with x = t/M 2 gives (/? = \ for the Bethe lattice) 



-21 a+ - + 4 J +o(t) 



(44) 



wnere a .- 4 (g _ 1)2 a c - 12c . 

The whole effect of having a = for the Bethe lattice can thus be absorbed into a 
redefinition of the constant a in the scaling function ((39) and (41)). Equation (36) can 
now be applied, with 

F(x) = d + b/x + c/x 2 (45) 
to obtain the scaling form of the Gaussian curvature on the Bethe lattice as 



b 2 1 fl 

8dt 2 (b + 6c/x) 2+ °\t 



2ql(q-2){K c t+<-K^AP}'' \t 
in agreement with equation (26). 
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(46) 



5 Conclusions 



In this paper we have explored the geometrical concept of a metric on the space 
of states of a thermodynamic system, two states are defined as being far apart if the 
probability of a fluctuation between them is small. In particular a class of models which 
exhibits a phase transition, Ising models on a Bethe lattice, has been examined. 

Using expectation values to define a metric on the space of states, as in equation (6), 
the geometry of the Ising model on a Bethe lattice with co-ordination number q has been 
analysed and the Gaussian curvature is given by equation (23), which is the first main 
result of this paper. This reduces to the previously known result for the one dimensional 
Ising model [7] when the co-ordination number q = 2. The curvature is seen to be positive 
definite for all temperatures and magnetisations and diverges at the critical point, but is 
finite at all other points (in particular at all other points along the curve of spontaneous 
magnetisation). The Gaussian curvature tends to the constant value q/2 for T — > oo as 
well as for M = ±1 (a particular case of the latter situation is, of course, T — > — for 
T > it can only be obtained in the limit of the external field h going to infinity). 

For a general ferro-magnetic the scaling form of the Gaussian curvature near the 
critical point has been calculated in terms of the scaling function and it's derivatives, (36), 
which generalises the t > 0, M = result of [2]. This is our second main result. It has 
been shown how this relates to the scaling form of the Gaussian curvature on the Bethe 
lattice, equation (23). 

An interesting open problem associated with the above analysis is to ascertain whether 
or not there is a relationship between the renormalisation flows for the Ising model on a 
Bethe lattice and the geometry described here. For the special case of the one dimensional 
Ising model it is already known that the renormalisation flow from T = oo to T = 0, along 
the M = axis, is a geodesic and that this is the only renormalisation trajectory that is a 
geodesic [6]. For the more general case with q > 2 it is certainly still true that the M = 
axis is a geodesic — this follows from the observation that the metric (and curvature) are 
even functions of M and thus invariant under a change in sign of M — and this axis will 
be a line of renormalisation flow, but it is an open question as to whether or not any other 
renormalisation trajectories are geodesies for the q > case. This question merits further 
investigation. 

The author wishes to acknowledge support from the Alexander von Humboldt founda- 
tion as well as the hospitality of the Physikalisches Institut in Bonn where the manuscript 
was completed. Some sponsorship was also received from Baker Consultants Ltd., Ireland, 
networking specialists ( jittp: / / www.baker.ief) . 
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Fig. 1: Gaussian curvature for the Ising model on a Bethe lattice (co-ordination number 
q = 3). The curvature diverges at the critical point K = (1/2) ln(3) ~ 0.5493, M = 0, but 
is finite everywhere else in the stable region of the K — M plane. The line of spontaneous 
magnetisation corresponds to the left hand edge of the graph (the vertical axis has been 
truncated at 1Z = 100). 




o o o o o o , 

O 03 ID CN 

i — I 

Pi 



12 



Fig. 2: Comparison of the line of spontaneous magnetisation and the pseudo-spinodal line 
for the Ising model on a Bethe lattice with co-ordination number q = 3. The vertical axis 
is the magnetisation M and the horizontal axis is z = e~ 2K . The leftmost thick line is 
the pseudo-spinodal curve, along which the susceptibility diverges (this lies in the unstable 
region), the rightmost thin line is the line of spontaneous magnetisation, which corresponds 
to the left hand edge of the graph in figure 1. These two lines kiss at the critical point 
M = 0, z = 1/3. 
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